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96 LING. ON THE SOLUTION OF A CERTAIN DIFFERENTIAL EQUATION 

I. 

Introductory. 

1. Objects of tke paper. In his discussion of the kinetic theory of tides, 
Laplace found that the function expressing the height of the tide at a given 
point due to the attraction of the disturbing body satisfied a certain differen- 
tial equation. Finding himself unable to obtain the general solution of the 
differential equation, he applied himself to the discussion of several particular 
cases which arise when certain assumptions are made regarding the physical 
constitution of the ocean. One of the cases he treated was that of the semi- 
diurnal tide when the depth of the ocean is supposed to be constant. In the 
course of his treatment of this case certain considerations enter which have 
given rise to much discussion. It is proposed to devote some attention to this 
case, and it is hoped to extend the treatment of this case so as to include ■some 
phases of it not previously treated. While it is generally conceded that the 
facts in regard to the disputed point referred to, have been made evident, yet 
the methods of placing those facts in evidence have been called into question 
by several writers on the subject, and do not appear to be the most satisfac- 
tory ones that are available. 

II. 

Historical Sketch. 

2. Origin of the problem. As just mentioned, the subject to be discussed 
was first treated by Laplace. His kinetic theory of tides is set forth in the 
M^canique Celeste, and the part with which we are concerned is to be found 
in Livre IV of that work, his solution of the differential equation being given 
in Article 10. Considerable time had elapsed between his first discussion of 
the subject and the publication of his great work. The earliest presentation 
of his treatment of the subject was contained in a memoir* presented to the 
Acad^mie des Sciences, and contained in Tome IX of the (Euvres de Laplace. 
He has sought a solution of the equation in the form of a series of positive 
entire powers, and has made use of a certain infinite continued fraction in the 
evaluation of one of the coefficients of the series. The correctness of the value 
found by his method has been questioned. As the solution in the series form 
was made the basis of his calculations, it was of great importance that no mis- 
take should be made in the determination of the coefficients, and more espe- 
cially in the determination of those occurring early in the series. 

* Recherches sur plusieurs points du Systeme du Monde. Memoires de l'Academie royale de 
Paris, annee 1775-6. 
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3. Previous contributions to the subject. In his early memoir Laplace 
has gone somewhat more into detail, and the method by which he determined 
the value of the coefficient is clearly shown. In the later work he has omitted 
a great part of the explanation, and has contented himself with expressing the 
quantity in the form of the continued fraction to which reference has been 
made. The later presentation of the subject has been the more accessible of 
the two, and on it all the later writers appear to have based their remarks con- 
cerning Laplace's method, while the original presentation has been overlooked. 
Attention has been called to it by Prof. Lamb in his recent work on Hydro- 
dynamics, and to him seems to be due the rediscovery, so to speak, of the 
memoir. Laplace's evaluation of the coefficient was objected to by Sir G. B. 
Airy,* and later the same objection was made by Mr. William Ferrel. t A 
defence of Laplace was made by Lord Kelvin in the Philosophical Magazine 
for September, 1875. The October number of the same journal for 1875 con- 
tains a note written by Airy in which he reaffirms his objections to Laplace's 
result, and appears not to regard Kelvin's reasoning as convincing. The num- 
ber of this journal for March, 1876, contains a reply by Ferrel to the arguments 
of Lord Kelvin. Prof. G. H. Darwin in the Encyclopedia Britannica| gives in 
more detail Lord Kelvin's argument. His treatment of the subject may also 
be found in Basset's Hydrodynamics, Vol. II, and Basset has briefly referred 
to the subject in a foot note.§ The latest contributions to the subject are 
believed to be the two papers by Ferrel which appear in Volumes 9 and 10 of 
Gould's Astronomical Journal. The latter of the two papers may also be found 
in the collection of papers || on the " Mechanics of the Atmosphere " edited by 
Prof. Cleveland Abbe. Reference may also be made to Professor Lamb's 
Hydrodynamics, in which attention is called to Laplace's original memoir. 

4. Synopsis of previous contributions. Before treating the problem 
analytically it will be useful to sketch the arguments of Laplace and those who 
afterwards treated the subject. Laplace, assuming that the solution of the 
equation could be expressed by means of a Taylor's series, substituted such a 
series with undetermined coefficients in the differential equation, and was able 
to determine all the coefficients of the series in terms of one of them, which 
remained arbitrary. He had previously argued that it was not necessary to 
obtain the general solution of the equation, since, as he affirmed, the arbitrary 
constants would be determined by the initial conditions of the water and would 
introduce effects dependent on this initial condition, which effects ought to be 

* Article, " Tides and Waves," Encyclopedia Metropolitana. 

t " Tidal Eesearches," Appendix to United States Coast and Geodetic Survey Eeport for 1874. 

X Article, " Tides," Encyclopedia Britannica. 

§ Basset, Vol. II, p. 218. 

|| No. 843, Smithsonian Miscellaneous Contributions. 
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disregarded, since in the case of the sea they would long ago have been over- 
come by friction. Considering, then, that any particular integral was suffi- 
cient, he proceeded to choose the most satisfactory value of the coefficient. 
His method of deciding the proper value of the coefficient will be given in 
Section III. It enabled him to satisfy himself that a comparatively few terms 
of his series would give the result with a very small error. From the form in 
which the result is set forth in the M^canique Celeste it appears, however, 
that he made the assumption that the ratio of any coefficient in his series to 
the preceding one becomes ultimately smaller than any assignable quantity. 
Moreover Laplace's argument regarding the sufficiency of any particular solu- 
tion did not occur immediately in connection with the treatment of this par- 
ticular case, and it therefore appeared that he offered no justification for the 
assumption. Airy objected to the assumption on the ground that it was 
unnecessary and unduly specialized the solution. He added that, if the sea 
were bounded by a parallel of latitude instead of covering the whole earth, 
then the arbitrary constant could be determined from the corresponding 
boundary condition. Ferrel agreed with Airy, and regarding the constant as 
being entirely at his disposal, based his calculations on the series resulting 
from assigning to it the value zero. Kelvin in his reply to these arguments 
quoted Laplace's reasoning regarding the sufficiency of a particular solution, 
but pointed out that this reasoning was not correct. Proceeding, he contended 
that, as demanded by Airy, there was a certain physical condition to be satis- 
fied, and that this condition was sufficient to justify Laplace's result. He also 
pointed out that the general solution of the equation should contain two arbi- 
trary constants, and that a further boundary condition would be necessary for 
the determination of the second of these constants. He showed that, since the 
oscillations of the water which are taken account of in the differential equation 
have a perfectly definite period depending on the period of the disturbing 
body, the original state of motion could not be taken account of in the solu- 
tion, for, except for special depths of water, the period of the latter oscilla- 
tions would be different from that of the former. Airy and Ferrel, however, 
did not admit the force of the reasoning by means of which Kelvin justified 
Laplace's result, and Ferrel's later papers are devoted to an attempt to show 
that the determination of the value of the constant is unnecessary. While it 
would seem that the constant was correctly determined by Laplace, it appears 
to the author that the analytical proof of this fact indicated by Lord Kelvin, 
and given in greater detail by Prof. Darwin, is not complete. It seems, too, to 
be desirable to obtain the general solution of the differential equation, and to 
follow up the suggestions of Lord Kelvin regarding the application of this 
solution to the more general case and some special cases. 
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III. 

Laplace's Solution. 

5. General outline. Apart from physical considerations, the arguments 
made by Laplace from analysis do not appear to offer a very good reason for 
his evaluation of the coefficient. In order to show this, it will be necessary 
to give Laplace's discussion as it appeared in his original memoir. The form 
of the equation as treated by Darwin differs slightly from that in which Laplace 
used it, but no essential difference is introduced by the change, and the same 
difficulties arise in both cases. As Darwin's form of the equation is doubtless 
that in which the equation will henceforth be studied, it has been adopted 
here. The equation may then be written 

a? (1 - a?) U - x ^ - (8 - W - fr) u + £fh* = 0. (1) 

It is to be noted that x is the sine of the polar distance of a particle of water, 
u the difference between the tide height in the dynamical theory and the tide 
height in the equilibrium theory. Assuming as the solution of (1) a Taylor's 
series containing only even powers and with undetermined coefficients, Laplace 
found that the coefficient of x* remained undetermined. He next proceeded 
as follows : He assumed as an integral a sum of a finite number of powers, 
and found that by adding a certain term to the left member of (1) this equation 
could be modified so as to have as a solution the assumed function. By study- 
ing the effect of increasing the number of terms in this function, he came to 
the conclusion that a very small error would be made in assuming as a solu- 
tion, such a function with a large number of terms. 

6. Application. To apply this treatment to the equation (1), assume 

' Rt « = A + (A, - E) * 2 + S k Aj?« . (2) 

2 

If this function satisfies equation (1), the following relations must be satisfied 
by the coefficients : 

A = , (a) 

A x -E=0, (b) 

A 2 =A it (c) 

16^3 — 10A 2 + jiA l = 0, (d) 

2A k+1 [2 (A - If + 6 (A - 1)] - 2 ^[2 (A - If + 3(A - 1)] + /Li*-i = , (e) 

(A = 3, 4, 5, ... , r.) 

- 2A +1 (2r 2 + Br) + (3A r = , (f) 

+ /3A+i = . (g) 
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There are, since (c) is an identity, r -\- 3 linear equations to be satisfied by 
r -\- 2 unknown quantities. It is easy to see that all cannot be satisfied ; for, 
starting from (g) and working back, there result zero values for all the A's and 
this disagrees with (b). If, however, one of the equations be rejected, the 
remaining r -f- 2 relations are sufficient to determine the values of the A's. 
Suppose (g) to be rejected. Following Laplace's method let the following 
abbreviations be made : 

,x r _, = 2(r - l) 2 + 3(r - 1) - [(r - l) 2 + 3(r - 1)] /i/ft , 
p r _ k = 2 (r - £) 2 + 3 (r - A) - [(r - A) 2 + 3 (r - *)]>/^*+i , 
Then it follows that 



whence 



-^r+l — 


ft "' 


A = 


ft— 1 


A-l = 




A = 


ft ft 


A 3 = 


ftft 



ftftft, 



A+l= ^ ^. 

ftftft • • • ft 

These values of J. 1( -4 2 > ^3> • • • > A r+1 satisfy the equations (a), (b), (c), (d), 
(e) (f), but equation (g) is not satisfied, and the function (2) is therefore not a 
solution of (1). If, however, to the left member of (1) be added the quantity 
— plA r+1 x i '' +6 , equation (g) becomes the identity 

- {3A,. +1 + pA r+1 = , 
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so that the expression (2) is a solution of equation (1) thus modified. The 
new equation can be written 

a? (i _ a?) |J _ x p - u (8 — 2a? — faf) + Epx* ?£— «?'+« = 0. (3) 

Laplace then argued that if the corrective term were very small, only a 
small error would be made if (1) were replaced by (3). He then proceeded to 
show that by taking r large enough the corrective term could be made to 
decrease indefinitely. 

7. Objections to method. In order that the discussion just given may 
justify the choice of the value of A 2 , it should put in evidence some property 
possessed by the series when Laplace's value is given to A lt and not possessed 
by it under other circumstances. All that is attempted in the preceding process 
is to show that the error made in assuming as the integral a finite number of 
terms belonging to the infinite series can be made less than any assignable 
quantity. But this is true of any series which is an integral of the equation 
and whose region of convergence is large enough to suit the conditions of the 
problem, and it will appear that no matter what value be given to A % the region 
of convergence of the resulting series is still large enough to suit the purpose. 
Moreover, it is not definitely proven that a vanishing correction to the equa- 
tion necessarily indicates a corresponding vanishing of the error due to the 
assumption of a finite number of terms instead of an infinite number. 

8. Relation between correction to equation and error in series. The rela- 
tion between the correction to the equation and the error in the value of the 
dependent variable can be shown perhaps more clearly in the following man- 
ner, assuming certain properties of the series used which will be deduced later : 
The differential equation (1) can be regarded as a linear relation connecting 

(7 it d it * 

the quantities u, -^- , and -j-^ in which the coefficients are rational entire func- 
CvOj ecus 

tions of a?. The function u is to be expressed by means of an infinite series. 

This series will have a certain circle of convergence. For all points within 

fffi ft U 

this circle the corresponding series for -j- and -t-j will also converge. Sup- 
pose the circumference of this circle of convergence lies entirely outside of the 
boundary of the region in which the independent variable is to vary. Then 

when w, ^ , and -^ are each replaced by the finite number of terms from the 

series expressing their values, certain errors will be made in the case of each. 
These errors will each become less than any assignable quantity if the number 
of terms be sufficiently increased. Let (1) then be written 

d 2 u . a du . » 

a lM + &dx + rU = d > 
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where a, /?„ y, and d are finite for all values to be considered. Let the true 

(J 1J ft U 

values of -j-^ , -=- , and u when found from the infinite series be A, B, and 

O ; and let s„ s 2 , and e 3 be the corresponding errors made in taking the finite 
number of terms for each of the three quantities. Then 

aA + &B + yC = d , (h) 

and 

a (A + e,) + &(£+ e 2 ) + y{C + t s ) > 3. (k) 

Then the correction to be added to the left member of (k) to make it an 
equality is — (ae, -j- ^e 2 -f- ye 3 ). Under the circumstances assumed above, in 
reference to the region of convergence, this correction will become indefinitely 
small when the number of the terms is sufficiently increased. But other cases 
may occur. Consider the case of the infinite series in which A 2 is given auy 
other value than that assigned by Laplace. It will afterwards appear that this 
series converges for all values of x within the unit circle, and also for * = ± 1. 

ft ft (it J 

The series for -j- and -j-g , however, are convergent only for points within the 

unit circle. It is clear then that, for all values of x less than unity, e u e 2 , and 
e 3 can be made less than any assignable quantity, and that therefore the same 
is true for ae lt -f- ^s 2 , -\- ye s . 

But x must be considered for all values up to and including unity, and it 
does not appear that, as x approaches unity, the correction must necessarily 
indefinitely diminish, since e 3 is the only one of the three quantities which 
indefinitely diminishes. Moreover the series is a satisfactory solution if only 
e s can be made indefinitely small for x = 1. On the other hand, it would also 
appear that the correction might be evanescent when e 3 is not so. 

9. Assumption equivalent to Laiylacds assumption. The same set of 
equations for the determination of the coefficients will be obtained if it be 
assumed that 



£ A r+1 
A r 



ft 1 = 0. 

A r 

For then it would be proper to assume that the equations connecting the coef- 
ficients eventually took the same form as (f). For all other series in which the 
relation just written is not true the equation (f) is not satisfied when a finite 
number of terms is taken, but has to be corrected by the addition of a term. 
From this point of view then Laplace's process would seem to necessarily lead 
to a series convergent over the entire plane. 
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IV. 

The Solution of the Differential Equation. 

10. Character of the integral. Laplace's solution having been considered, 
the general integral of the equation, 

* {1 -^ S - X Wx ~ U{8 ~ 2 * 2 _ ^ = -P^' (1) 

may now be sought. In this equation it must be remembered that x = sin 
where d is a polar distance. It is necessary first to solve the auxiliary equa- 
tion, 

x*(l-x>)^-x^-u(8-2x?-px*) = 0. (7) 

The following general theorems will be useful : 
Theorem I. In order that the equation 

d n u . m d n ' l u . m d n ~ 2 u . d n ~ 3 u du , A ._, 

"3SF +Pi d^ + ^d^ + ^dx^ + --- + *^dx +P« U =° W 

shall have n independent integrals of the form, 
u=af (P n _ x log-"'* + P„_ 2 log"- 2 * + P n _ 3 log"" 3 * + ■ ■■ + P 1 \ogx + P,), 

where P , P lt P 2 , P 3 , . . . , P„_ x are expressible in the neighborhood of x = 0, 
in series of positive and negative powers of x, the number of negative powers 
in each being finite, it is necessary and sufficient that for each of the coefficients 
of the equation, such as p it the point x = shall be an ordinary point, or a 
pole whose order of multiplicity does not exceed i. 

These integrals will constitute one or more groups of the form 

u x = x r M x , 

u 2 =z x r (M 2 log x + iig , 

u 3 = x r (M 3 log 2 * + JV 3 log x + i? 3 ) . 



u k = x*- (M k log*" 1 * + N„ log*- 2 * +...), 

where M x , M 2 , M 3 , . . . , M k differ only by constant factors. 

Theorem II. The integral of the equation (8) will be continuous and 
monogenic for all values of x for which the coefficients p u p 2 , p 3 , . . . , p n are 
continuous and monogenic, and it can possess no critical points which are not 
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also critical points of one or more of the coefficients. It may not have critical 

points at all the critical points of the coefficients. 

The proofs of these theorems may be found in Jordan's Cours oV Analyse.* 
If, then, equation (7) be put in the form (8) and these theorems applied, 

it is clear that, 

(1) Equation (7) has one integral which can be expressed in the form of 
a power series <[>, and a second integral which can be expressed in the form 

^2 + <P\ l°g x > where (p % and <j> x are expressible as power series and ^ = a con- 
stant c (when c = 0, the second integral is also expressible as a power series) ; 

(2) The general integral of (7) can have no critical points except at 0, 
± 1, ± oo . 

(3) If x be replaced by — x, the equation remains unaltered ; whence it 
follows that the function has the same character at x = — a as at x = + a, 
and, in particular, if the general integral of (7) has a critical point of any sort 
at x = + 1 it has a critical point of exactly the same character at x = — 1. 

It follows, too, that if U is any particular integral of (1), the complete 
integral is 

u = a x <}> + a 2 (</> 2 + c»p log x) + U, 
or 

u = a x (p + a 2 </> 2 + U 
when c = 0. 

11. Deduction of the complementary function. It is correct, then, to 
assume as one integral, 

u = I,. A,x m+r ° , 



s being a positive integer. Then 

^' = 2; (m + rs) A r x m " 1+rs , 
ax o 

~ = 2 r {m + rs) (to — 1 + rs) A,x m ' 2 + rs , 

The result of substituting in (7) is 

2V (m + rs) (to — 1 + rs) A r x m+rs — I r (to -f rs) (to — 1 + rs) A r x m+2 +'' a 



— I,, (to + rs) AtX™*'' — 2 r 8A r x m+rs 



+ X 2A r x m+2+r * + X 0A r x m +*+ r ' = 0, (9) 



* Vol. Ill, Arts. 146, 92, 118. 
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identically, The coefficient of x m , the lowest power of x, must vanish ; whence, 
since A o 0, 

?n? — 2to — 8 = 0, and to = 4 , — 2 . 

The value of s must be 2, for if s o 2 the coefficient of jc" 14 " 2 does not vanish. 
The result of substituting s = 2, to = 4 in (9) is 

l r (4 + 2r) (3 + 2r) ^ r+4 — J r (4 + 2r) (3 + 2r) ^ r a; 2r + 8 



— f r (4 + 2r) ^l r a.>+ 4 - 1,. 84 r * 2r + 4 



+ l r 2^1 r « 2, -+ 6 + l r S^^+s = . (10) 

' 

From (10) the equations for the determination of the coefficients are, 

16^, - 10A = , (11) 

2k (2k + 6)A h ~ 2k (2k + 3) A k _, + (3A h _ 2 = , (12) 

(* = 2, 3, 4, 5, . . .) . 
Thus each coefficient is a multiple of J. and one integral can be written 

«, = a? 4 + 4 1 * 6 + 4*^ + 4- 3 * 10 + • • • 

A„ jl _4 

(13) 

= x* + Off + C#? + Ctx™ + .... 

It is easy to verify that K± u x is the series written down by Airy in his 
article in the Philosophical Magazine as the correction to Laplace's series, K K 
being an arbitrary constant. On applying to the equation the theory set forth 
by Heffter,* it appears that there is no series corresponding to the root 
to = — 2. The second integral is then obtained in the form 

U 2 = 02 + ^1 l g X = <Pz + -A.' "i log * . 

From this 

dx ~ dx +JL l0gX ^ + l^ Ul ' 

dx*~ dx* + x dx F ' + g W • 

* Lineare Differential gleiehungen, § 16. 
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These values being substituted in (7), there results 
^(1-^)^-^11-^(8-2^-^) 

— A% (2 — a 2 ) + A'.2x (1—x 2 )^ 



+ A' log x 



* 2 (1 ~ J)^?-"^ ""I ( 8 - ^ ~ M 



, (15) 



identically. If the substitution (p 2 — A'ip 3 be made and it be borne in mind 
that Mi is a solution of (7), (15) reduces to 

*0--<*)%-*^Z ~ &(8 - 2*" - /to*) 

-(2-x 2 )v l + 2x(l-x>)^ = 0. (16) 

If now it be assumed that 

4> i = S r B 1 ^+", (17) 



and the substitution be made in (16), the equation for the determination of m 
is, as before, 

m 2 — 2m — 8 = 0, whence m = 4, — 2 . 

The theory shows that m = — 2 is the root to be taken, and that, as before, 
s = 2. The substitution in (16) gives 

1; (2r — 2) (2r — 3) B r x 2r ~ 2 — l' r (2r — 2) (2r — 3) B r x* 



— I r (2r — 2) BrX*- 2 — 1, 8B r a? r - 2 + I r 2B r x 2r 



+ 1,. pB r x»-+ 2 — i',.2(7 r « 2 '-+ 4 + 1; C r as*+« 



+ 1',. 4 (r + 2) C^ 4 — 1; 4 (r + 2) tf r a; 2r + 8 = , (18) 



in which the (7's are the coefficients in the series u t . 
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From (19) are derived the relations 

= 8B X + 4B , 

= 8B 2 -2B x -pB , 

6= -fa, 

(19) 
7 - IOC, = 16^ — 10^ + P£* , 

(4k — 5) C k _ 3 - (4k — 2) <7*_ 2 = (2k — 4) (2k + 2) B k+X - (2k - 4) (2k - l)B k 

+ fiBk-i • 
(k = 4, 5, 6, . . .) 

These equations determine all the coefficients in terms of B s , the coefficient of 
a; 4 , which remains arbitrary. 

Since u x is a solution of (7), it is clear that, if <p is a solution of (16), so 
also is f + A' 3 u x , A' 3 being any constant. Moreover u x starts with the fourth 
power of x. Then if any value be assigned to B 3 and the resulting value of 
<p 3 be denoted by f we can write 

& = f + A\u x . 
<p is a series of ascending entire powers starting from B^e' 2 in which every 
coefficient is known. B 3 may, if it is desired, be taken to be zero. The com- 
plete complementary function of (17) then is 

u = A\ (fa + <p x log x) + A\u x 

= A \ (& + A ' u \ !og «) + -4>i 
= A'A\ (<j> 3 + u x log x) + A\u x 
= A (<p + u x log x) + Bu x . (20) 

12. The particular integral. It remains to determine a particular inte- 
gral of the complete equation (1). The character of u x and of the absolute 
term of (1) makes clear the existence of a particular integral expressible in the 
form of a series of positive entire powers. Assume then 

U = A B + (A x — E) x 1 + I A^ . 

After substitution in (1) there result the relations 

A = o, 

A x — E=0, 
A 2 = A 2 , 
2A k+x {2(k— If + 6"(A — 1)} -2A k {2(k-Vf + 8 (A - 1)} + £>A k _ x = 0. 

(k = 2, 3, 4, . . .) 
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All the coefficients are given in terms of E and A 2 , the latter being arbitrary. 
Since only a particular integral is required any value may be given to A 2 . In 
proceeding to a choice of a value it is interesting to follow the method by 
which Laplace's continued fraction is obtained. From the relations written 
just above, it can easily be deduced that 



-4*+i 



: 2 (2k* + 3k) - 2 ('2k* + 6*) A k+2 /A k+l " {k ~ % 3 ' 4 ' ' ' " } 



A 2 _ p_ (2.1 2 + 6.1)/3 (2.2 2 + 6.2)/3 

A l 2(2.1 2 + 3.1)— 2 (2. 2 2 + 3.2) — 2(2.8*"+ 8.8)"— '" ' 



{2(« -1) 2 + 6( »- l)}j9 



2 [2 (» - l) 2 + 8 (n - 1)1 — 2 [(2w 2 + 8») - (2n? + 6n)] yl^/A 



and 4, = i?. 

The assumption, 



£ 



71=00 



4^ = ° . (21) 



gives the value of A 2 in the form of an infinite continued fraction. It is per- 
missible, if convenient, since any value of A 2 may be taken. It was made by 
Laplace in the M^canique Celeste apparently without justification ; but, as 
has been seen, Laplace believed in the sufficiency of a particular solution and 
considered the resulting series as a satisfactory solution without the addition 
of the complementary function. The assumption (21) so affects the coeffi- 
cients that the series converges for all finite values of x. It is not necessary 
that the particular integral should converge for points outside the unit circle. 
It is convenient from a mathematical point of view, however, to choose this 
series as the particular integral, for, if it were necessary to study the function 
for points outside the unit circle, it would be sufficient to obtain the comple- 
mentary function in the form of a Laurent's series while the series just found 
would serve again as a particular integral. The assumption (21) is seen to be 
equivalent to that involved in Laplace's original process. From the reasoning 
of this section it is clear, too, that when 



then 



£ -^n+l — > 
~A\~< ' 



£■'4-1+1 __ j 
A n 
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From this point of view, too, it is clear that the series under consideration 
converges at least for every point within the unit circle, and that if it converges 
for a greater circle of convergence it converges for every finite value of x. 
For the sake of definiteness Laplace's value of A 2 will be denoted by L and 
the series which furnishes the particular integral of (1) will be denoted by 2. 

13. Properties of the complete integral. The complete integral of the 
equation (1) can then be expressed for points within the domain of the origin 

by 

u = A (<p + u x log x) + Bu x + 2 . (22) 

(1) Convergence : 

The most general integral in the form of a positive power series can be 

written 

u = B u x -f- 2 . 

The relations among the coefficients of such a series show that, unless B 
is zero, the circle of convergence is of unit radius ; and when B is zero, the 
circle of convergence has an indefinitely great radius.* It follows, then, that 
u x converges only for points within or on the circle of unit radius. Again 
(p 3 = f + A 2 u x does not converge for points outside the unit circle. Then <p 
converges for points within the unit circle. It is conceivable that B 3 may have 
been chosen so that f shall converge all over the finite part of the plane. 

(2) If A = B = 0, the function has no critical point except at infinity. 
If A = 0, B < 0, the function has critical points at ± 1, ± oo . If B = 0, 
A ^ 0, the function has a critical point at 0, ± oo , and (except for one par- 
ticular choice of B 3 ) at ± 1. In addition to the singularity of <p at x = 
integrals of this class and integrals of the general class have a singularity at 
x = due to the singularity of log x at that point, and are, in addition, many 
valued at any point owing to the properties of log x. This indeterminateness 
will be removed if it is assumed that for positive real values of x the result 
shall be real. Again, when B =: 0, 



du - 

dx 

But du du 



'£+**•£ 



,1 T , d2 



dd ~ dx vl ~ ** : 



£fdu~} f» (du .- -„ 

{do\^rl{^ vl - 3? 

3=1 

= a£[ </T=l# d £] + a£ [(log x) VT=r^ 



du x 

dx 



' The detailed proof of this faot is quoted in Section VI. 
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It will afterwards be shown that for the complete integral (22) 

= a finite quantity, 



and that 



du 

To 



Cdu x ' 



= a finite quantity. 

9=ir/2 



.•. for all integrals of the form 

u = A (f -\- u x log x) + £ , 

I -t^M = a finite quantity or zero. 

V. 

The Determination of the Constants for Laplace's Case. 

14. The physical conditions. It having been agreed that the constants 
shall be determined to suit the boundary conditions, the case discussed by 
Laplace, where the whole earth is covered with water, may now be treated. 
The expression for u, as given in (22), has an infinite value when x = 0, unless 
A = 0. Since there cannot be a tide of infinite depth at the pole it is neces- 
sary to make A = 0. The remaining expression is 

n = Bu v -f 8 . 

Airy and Ferrel contended that this was the exact expression for u, and that 
B could be given any value. Ferrel determined it by the condition 

B + L = . (23) 

Lord Kelvin pointed out that owing to the symmetry of the disturbance in the 
two hemispheres the meridional displacement of water should vanish at the 
equator. The expression for the meridional displacement is the product of 
two terms of which one involves the latitude and the other does not. The 
factor involving the latitude is, 

b -- -.~~ ,r, % + 2m cot d] . (24) 

du j 



4m sin 2 
When = jt/2 * 



_ J^ (du 

Z ~ 4m [d0_ 



«-ir/2 



* = co-latitude or polar distance. 
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at the equator it is necessary that 

du 
dd 



. (25) 



But 



du .- -^ du 

dd~ Vl ~ af dlc' 



it is necessary that, 



£{ 



^ TZ ^£\=«- w 



15. Proof that B = 0. In order to prove that the condition (26) requires 
that B shall be zero, the function will be considered in the neighborhood of 
x = l. 

Assume 

x --■= 1 + y ; (27) 

then 

du du d 2 u d 2 u # 

dx dy ' dx* dy' 2 ' 

and equation (7) takes the form 

(2y + 5tf + if + y*) + (1 + y) g 

+ «{(6 -ft - 4(1+ fiy -2(1 + mf-±Pf-W} =0. (28) 
Theorems 1 and II apply to (28) also. Then assume 

u = l' P A r y m+r > , 



s being a positive integer. The equation for the determination of m is 

2m 2 — m = ; (29) 

whence m — or + J. Also a = 1. The two integrals are expressible in 
series form. Moreover the relations connecting the A's are, for both values 
of to, such that each coefficient is given as a multiple of the first one. The 
two independent integrals of (28) may then be written 

y x = 1 + 2' r a r y r , 

(30) 

i 
the «'s and /?'s being known. 
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The complete integral of (28), then, is 

M = Wi + c 2 y 1 . (31) 

If now the substitution (27) be made in (1), an equation results which differs 
from (28) only by the expression — Eft (1 -\- yf in the right hand member. 
Of this equation, the complementary function is given by (31). To complete 
its integration, it is necessary to find a particular integral. The character of 
y v and of the absolute term — Eft (1 -f- yf make it clear that a particular 
integral can be obtained in the form of a series of positive entire powers. 
Let this integral be denoted by Y It is not of importance that its coefficients 
be calculated. Then, in the neighborhood of x = 1, or of y = 0, the integral 
of (28) can be expressed in the form 

u = C\9i + c& 2 + Y. 



Then 
Also 



dx dy l dy dy dy 



l/l - a? = V- {f + 2y) = iyii/2 + y , 
where i = j/ — 1 . 

••• m = d £ ^r^ 2 - <w» (2 + yf f + tf * (2 + yf f 

+ ¥(2 + 2/)»^. (32) 

When y = 0, the first and third terms on the right vanish, and the middle 
term becomes 

[+ * ST* cj,yi (2 + y)*]^, or c 2 i/V2 ■ 

(du/dd) e=n/i = c 2 i/}/2 . (33) 

In order, then, that at the equator du/dO = 0, the function must be such 
that c. z = 0. But if c 2 = 0, the point x = 1 is an ordinary point of the func- 
tion, while if c 2 o the point x = 1 is a branch point of the function. If the 
point x = 1 is an ordinary point of the function, so also is x = — 1. 

In this case the function has no critical point in the finite part of the 
plane, and, if expressible in the neighborhood of the origin x = by means 
of a Taylor's series, that series will converge for all finite values of the varia- 
ble x. If the point x = 1 is a branch point of the function, and if the function 
is expressible in the neighborhood of the origin x = by means of a Taylor's 
series, that series will have a circle of convergence of unit radius. It follows, 
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then, when B = and u = 8, that c 2 = and {du/dd)%~ n /i = 0, and when 
-6^0, that c 2 ^ and (du/dd) t=n/i ^ 0. 

Thus it is seen that the condition stated by Lord Kelvin requires that 
B = and u = S. Consequently the series as written by Laplace is the 
complete solution for the case of an earth completely covered to a constant 
depth by water. 

Note. — Equation (33) gives an imaginary value for {du/dd) 6=v/i when the 
arbitrary constant c 2 is taken real. It is evident, however, that, when real 
values of the function between x — and x = 1 are desired, c 2 must be taken 
purely imaginary ; for it was assumed that 

x = 1 + y ; (27) 

so that when x is less than unity y is negative, and y* is a pure imaginary. 
The cg/l will be real when c 2 is purely imaginary, and it follows that then y 2 is 
also real. 

VI. 

Darwin's Presentation of Lord Kelvin's Proof that B most be Zero when 

(du/dO) e=n/i =0. 

16. Darwin's argument. The function u = Bu y -f £ may be regarded 
as a single series of even positive integral powers commencing with the fourth 
and having the coefficient of x* arbitrary. It has already been seen what rela- 
tions connect the coefficients and define them in terms of the coefficient of x* 
(A 2 say). It is known, too, that when 



then A, = L. 



2 



Suppose now that 



£ 

£ 



^±j = o 



A n < U ' 



but = a -1 , a finite quantity. Then 



A n+2 2n + 3 
A n+1 2n + 6 


2n(2» + 6)A n+1 


2n + 3 
2n + 6 


2n(2n + 6)^ 



(34) 
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where h tends to zero when 'n becomes indefinitely great. Darwin's argument 
is along the following lines : When 



•^n+l > ft 



then, for large values of n, 



A„+ 2 X 



,2n+4 



A ~2n+2 



> + 3 a* _ 
2n + 6 



2 (n + 3) 



Z 2 



fi- 



L 



2w 



tf 2 



nearly. 

But if (1 



a; 2 )* be expanded by the binomial theorem the ratio of the 

3 1 



(n -f- l)th term to the nth term is 
it is possible to write 



1 — 



2n 



3?. Consequently, in this case, 



u = A l + B l (l — x*)i , 

where A l and B x are finite for all values of x. A similar argument being made 
in the case of the series for du/dx it follows that 

du/dx = C + D(1 — x*)-* , 

where C and D are finite (and not zero) for all values of x. 
But 

du/dd = du/dx (1 - a?)* = <7(1 — x*)l + D ; 

.: (du/dd) e=ny2 = [C(l - aty + Dl^ = D>0. 

17. Discussion of Darwin's proof . These results agree with each other, 
and with what has been proven in another way ; but this proof of the fact 
that B x and D are not zero nor infinite does not appear to be entirely satisfac- 
tory, and it is essential that this property of B t and D be made evident. The 
ratio A n+2 /A n+x becomes 1 — f n~ x when the square and higher powers of n~ l 
are neglected. If after a certain value of n quantities of the order n~ l be 
neglected, the ratio A n+2 /A n becomes unity. Then, following a line of argu- 
ment similar to that given, it would appear that 

u = A 2 + B 2 (1 - a?)' 1 (a) 

and, by a similar course of reasoning, that 

du/dx = C 2 + D 2 (l — x 2 )- 1 . 

These results do not agree and are incorrect, but they show in what respect 
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the previous reasoning is weak. For, suppose that the binomial expansion of 
(1 — a? 2 )* be written 



Then, if the infinite series can be written in the form 

u = A, + B x (1 - a»)» , 
where A x and B x are finite for all values of x, it follows that 



Now for every finite value of n (/? being positive), 

-^■r+l <- Tr- 

A r y, 

£ -A r +i __ { Tr+l . 
A r X* Y,. ' 



Al+1 < Tr+i . (b) 



Darwin has shown that 



but it is necessary also to show that 



i A 2 A 3 


A r A r+l 
A r _! A r 


<& r\'n" 


Tr Tr+l 
Tr-1 Tr 



(c) 



is finite. 

Both numerator and denominator are zero so that the value of the quotient 
requires investigation. 

The discussion of (a) and (c) shows that in (a) A 2 and B 2 do not have 
finite non-vanishing values for all values of x. 

VII. 

18. Cases to be treated. It remains, then, to examine the other cases 
included in the solution obtained. Airy pointed out that in the solution of 
the form 

u = Bu x (x) + 8 (x) , 

B could be determined so that the solution would be suitable for the case of 
a sea forming a spherical cap and extending from the pole to an arbitrary 
parallel of latitude. Lord Kelvin pointed out that, if the general solution were 
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at hand, the two constants could be determined so as to obtain a solution suit- 
able for a zonal sea lying between two parallels of latitude. The most inter- 
esting cases to be dealt with appear then to be the following : — 

1. Case of a sea covering the whole earth. This is the case already treated. 

2. Case of a sea extending from the pole to a given parallel of latitude. 

3. Case of a zonal sea bounded by two parallels of latitude on opposite 
sides of the equator and. equally distant from it. 

4. Case of a zonal sea bounded by any two parallels of latitude lying in 
one hemisphere. 

5. Case of a canal lying along a parallel of latitude. 

Case 2. 

19. Polar sea. If the sea extends only to a given parallel of latitude 
from the pole, it is necessary that the meridional component of the motion 
should vanish for the corresponding value of x, the sine of the polar distance 
of the boundary. 

Then, as in Case 1, A =0, and the condition just named gives for the 
boundary value of x 

~ + 2« cot = . 
at) 

Let, then, 0, be the colatitude of the southern boundary, and sin 0, "== a,. 
Since 

du du a j a > r\ 

-,-r. = -j- cos a , and cos #, -> (J , 
do dx <- 

dx sin a 
Bu( («,) + £'(«,) + | { Bu x («,) + 8 («,)} = ; 

JS - «,<(«,) + 2m, («,)" ( ° } 

This gives the expression for u at any point in the form 

K = , ( ,,_«±m MlW . (36) 

In particular, lor the southern boundary, 

w(ai) = ax . « feJ<W-g(^W ^ . (37) 
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(36) and (37) give the amounts to be added to the tide deduced from the equi- 
librium theory. 

The total tide at any point then is 

U=u(x) + Fa?, (38) 

and at the southern boundary 

U = u («!) + Ea? . (39) 

As has been stated, FerrePs calculations were made for the series in which B 
was given by the relation 

B + Z = 0. 

Then, from (35) it appears that the tides calculated by Ferrel would be those 
existing on a circumpolar sea bounded by a parallel of latitude $k — 0, where 
a = sin d and satisfies the equation 

z= a g fa ) + ag(a) 

<*«! (a) -f- 2wi(a) v ' 

Case 3. 

20. Sea extending equally on both sides of equator. Suppose the sea to 
extend equally on both sides of the equator, the boundaries being parallels of 
latitude. 

The condition that there shall be no motion of water along the meridian 
at any point of the northern boundary gives one relation connecting A and B ; 
but it is clear that the corresponding condition for the southern boundary gives 
exactly the same relation ; so that one of the constants appears to be arbitrary. 
The considerations which applied to Case 1 apply to this case. The symmetry 
of the motion requires that there be no meridional motion of the water at the 
equator. In this case, also, it is necessary that 

(du/d8), =w , a = . 

This gives a second condition by means of which the remaining arbitrary con- 
stant may be determined. 
From 

u = A {<p + m, log x) + Bu x + $ (22) 

it follows that 

%t = A(l- z 2 )l (<p ' + u( log x + ± i h ) + B(l - x*)i u' + (1 - a?)» £'. 
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Now 



£[(1-^)*S'J=0; 



= o 



since u x converges for x — 1 ;* 

£[(l-x*)lu l 'logx]=0, 



since log 1 



= and £ 



[ (1 — a? 2 )* m/ ] is finite. 



It is clear, then, that 



£[(i-* 2 )V] + ^£ 



*)*& +^t,[(l-a?)*Mj=0. 



It has been seen that 



-V [(1 — x 2 )' m/] is a finite quantity, say b ; 
•V [(1 — x 2 )* f'] is a finite quantity, say a . 



(In one particular case it is possible that a might be zero.) Then it follows 
that 

Aa + Eb = . (40) 

Returning now to the condition first stated, let t be the colatitude of the 
boundary, and let sin 0, = a v It is necessary that 



(du , 2m 1 n 



I 



dx 



x J, 



(41) 



since cos 0, ^ 0. The resulting relation between ^1 and 5 takes the form 
A 



2 12 

¥ ("0 + - 9 («i) + < («i) lo 8 «i + ~ M i K) + - M i ( a i) lo 8 a 



r M' 



/(«,) + | Ml ( fll ) ] + £'(«,) + | £K) = . (42) 



* See equation (b) Section VI. The expansion of (1 — x*)i converges for x = 1. 
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An equation for u is obtained by eliminating A and B from (22), (40), 

and (42). The total tide is 

u + Ex*. 

Case 4. 

21. Sea bounded by two parallels of latitude on the same side of equator. 
Suppose the sea to be bounded on the north and south by parallels of latitude 
and to lie entirely within one hemisphere. 

It is necessary that, at the northern and southern boundaries, 

du 2 « 

dx x 

Let the boundaries have colatitudes 0„ 6 2 (6 l < 6 2 ), and let sin 6^ = a„ sin# 2 
= a 2 . 

Then the equations for the determination of A and B are similar to (42) 
and are 

r 2 12 

A f'(«,) + - ip («,) + <(«,) log «, + - ^< I (« 1 ) + - «i («i) log «, 



+ £ 



<(«0 + i «, (a,)l + «'(«,) + | 8 («.) = , 
«i J a \ 



r o 1 2 T 

^ I ?'(«») + - P K) + <(«i) 1( >g «2 + - Ml(« 2 ) + - «l(«2) log «2 I 
[_ "2 "2 "2 J 



K43) 



As before, 



4 ^ r< (a 2 ) + i «, («,)] + «'(«,) + | 8 K) = . 

L "2 J M 2 



(22) 



J. 0(a) 4 m, (a) log a] + Bu^x) + S(a) — u — . 
Then, eliminating 4 and _Z?, the equation for u is 

<p(x) + m,(k) log a; , w,(a>) , S(as) 

9 12 2 2 

$»'(«,) + - tp(a l ) + <(«,) log «! + — «,(«!) + - Mi(«i)log «i, <(«,) + - «i(«i), £'(«i) + - S(«i) 

9 12 2 2 

y>'(«2) + - f (« 2 ) + <(« 2 ) ^g «2 + - W l( a 2) + - «l(" 2 ) ^g «2, <(« 2 ) + - M l(«2)» S 'K) 4 - S(«2) 



«2 

2 12 2 
p'(«i) + - ^(a,) 4 <(«,) log a, + - i/,(«i) + - «i(«i) log a„ Vfa) 4 - «i(«i) 

«! «1 «1 "l 

o 12 2 

y>'(« 2 ) 4 - pK) 4 <(« 2 )log«a 4 - Mj(« 2 ) 4 - "i(«j) log « 2 , m,'(«2) 4 -«i(« 2 ) 

a 2 a 2 o^ «j 



(44) 
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(44) gives the value of u. The complete tidal expression is 

« + Ea?. 
The value of u at the boundary whose colatitude is t is given by 



<P\ a i) + ~ V( a i) + M '(«i) log a, + - «,(«,) + 1 ufa) log «„ «/(<*,) + - Ml (a,) 
?>'(«s) + - ?(«.) + "'K) log «2 + - «i(«j) + - Wi(aj) log a 2 , <(«,) + - w,(« 2 ) 

"2 "2 tt 2 #9 



f < a l) 



1 



f'( a i) + - M iK) 



«i(«i) 

<(a,) 



"(«j) 
, S(«.) 



; .(44A) 



f'K) + - <P(<h) + <K) log^ + i «,(«,) + I «!(«,) log - 2 , <(a 2 ) + ^ Mi(a2)i gv^ + 2 g , > 

"2 '*1 a 2 a 2 a l <^2 #2 



A similar equation may be obtained for the evaluation of u (a 2 ). 

Case 5. 

22. Canal of width 2d lying along a parallel of latitude. Suppose the 
zonal sea to narrow down to a canal of width 2d having as its northern bound- 
ary the parallel of latitude \n — d. 

All three of the functions in the expression (22) are expressible in the 
neighborhood of = d t in series form. Let them be expressed in this man- 
ner and let d be taken so small that powers of it higher than the first may for 
purposes of calculation be neglected. 

The second of the equations (43) then becomes, after simplification by 
means of the first one, 

f 2 2 2 1 

A I f"(«l) + - f'(«l) — —2 ?(«l) + - <(«l) + —2 «l(«l) 



+ B 



2 2 

+ |>l"(«l) + - «l'(«l) — — 2«l(«l)] ^g «l 

<(«.) + 1 <(«o - 1 2 Ml (« t ) 1 + srw + 1 8'(«0 - A *te) 



+ d(lA +mB + n)V\ — «i 2 = , (45) 
where £, m, n may be easily found in terms of «j. 
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23. Tide at point distant d from boundary of canal. For any point 
within the canal distant d from the northern boundary x = a, -\- d \/l — a, 2 . 
This substitution being made and the product dd being neglected the result is 

<P («0 , «i(«i) » S(«i) 



2 1 



, <(«>) + | ".(«.) , 2'(«i) + | «(«i) 



?"(«,) + 1 <p '(«0 + | <(«,) + -^ «,(«.)> <(«,) + £ «/(«,), «"(«,) + I S'(«.) 



a, a, 

f(«i) 



+ rf l/l — a, 2 



, M,(a,) 



, «(«,) 



2 1 2 2 

f 'fa) + - ffa) + - M i(«i). «i'fa) + - M ifa)> 2'fa) + 77 S ( a ') 



£ — «i log a, 



m 



+ <V1-«, S 



f'fa) + — M ifa) 
"i 

1 2 



, «,'(«,) 



. S'K) 



, «,'(«,)+?■ «,(«,), «'(«,) + ?8(«0 



?"fa) + 1 ?-(«.) + ? «/fa)+ A «.(«,), <fa) + 1 <fa), 8"(«o + ?■ «'« 



= ?< 



2 1 

f'fa) + - ffa) + - «i(«i) 



, <fa) + - «ifa) 



3 2 2 3 

?"fa) + - fX<h) + - Ml'(«l) + — j «lfa), "l"fa) + - <fa) 
"1 «1 "1 «1 

2 1 2 

p'fa) + — pfa) + - Wi(«i), <fa) + - Wjfa) 



+ M p/1 <Zi 2 <2 



«1 

£ — m log <*! 



«i 



m 



ai 



(48) 



It is easy to show that 

2 4 4 

I — m log a, = f '"fa) + - ?>"fa) — —2 f'fa) + ^3 f fa) 



+ - «Tfa) + — 2 M l'fa) — p W lfa) » 



to = «r(«i) + - vw — —2 M 'fa) + — 3 W lfa) » 



«l' 



a/ 



» = «'"fa) + | «"fa) - ^ 2'(«i) + A *(«0 • 
"i «i «i 



WHICH PRESENTS ITSELF IN LAPLACE'S KINETIC THEORY OF TIDES. 



123 



Equation (48) can then be put in the form 

2 1 

f '(«i) + - <f(0 + - M i(«i) 



<(«i) + -«i(«i) 



3 2 2 3 
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«^ «i «j «i 
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a i a i I "i 1 

o -i i 12 i 
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2 1 
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o o o 

y"'(«l) + - M l"(«l) — —2 M l'(«l) + —3 "l(«l)> <"(«l)» 8 "'(«l) 

*1 



«1 " " "' «l" ■ ' - ■ «!' 

The value of u at the middle of. the canal is obtained by putting 8 = d, which 
simplifies (49) somewhat. 

24. Canal of negligible width. If now the canal be taken so narrow that 
the width may be neglected, (49) is still further simplified. The results should 
coincide with those obtained by considering the motion as in two dimensions. 
In this case the equation for u is 



2 1 
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f"(«0+-« 1 '(a t ) — iwi(aO, <(«i). S "(«i) 



(50) 



124 LING. ON THE SOLUTION OF A CERTAIN DIFFERENTIAL EQUATION 

As before, the total tide is 

U=u+Ea?. (51) 

VIII. 

Summary of Eesults. 

25. Summary of I— IV. For convenience of reference, and in order to 
render the results available to any who do not desire to follow through the 
processes of obtaining them, it has been thought desirable that they should 
be restated in a separate section which, along with the historical sketch in 
Section II, would give a complete account of the state of the problem. Sec- 
tion III is devoted to a discussion and criticism of the analytical process by 
means of which Laplace obtained his value for the arbitrary constant in his 
solution. Objection is taken to the process employed for two reasons. In the 
first place, the reasoning used has not been shown to be, and does not appear 
to be strictly accurate. In this connection it may be said that in the exami- 
nation of the apparent inaccuracies it has been thought sufficient to indicate 
the weaknesses of the method rather than go into a minute discussion of them. 
The modern advances in the theory of Differential Equations make it appear 
probable that matters of this character will be treated differently -in future. 
In Section IV the complete solution of the equation is found, the expressions 
iuvolved being infinite series, whose regions of convergence are large enough 
to make possible the treatment of all cases that can arise. The regions of 
convergence of the series together with certain important properties of the 
integrals can be predicted from the form of the equation. The integral found 
is more general than that previously deduced involving the two arbitrary con- 
stants. The series used by Laplace enters this integral as a part of it. In 
the derivation of the integral the method of Laplace as given in the M^canique 
(J&este is made clear. In the closing paragraphs of the section certain prop- 
erties of the integral important in the application of the physical conditions 
are deduced. 

26. Summary of V- VII. Section V deals with the applications of the 
physical conditions to the determination of the arbitrary constants. One of 
the constants is immediately determined. The determination of the other 
involves greater difficulties. The analysis on which previous evaluations have 
been based is rejected and replaced by a determination of the value which 
appears to be entirely satisfactory. The condition which determines this 
constant is the condition stated by Lord Kelvin. The result of this section 
appears to be a complete justification of Laplace's series. 
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In Section VI the objections taken to the proof previously given in the 
determination of the second arbitrary constant are set forth. 

The last Section contains the discussion of five important cases in the 
theory of tides. The arbitrary constants in the general integral of the differ- 
ential equation are determined so that the integrals represent the tidal disturb- 
ance in these cases, and expressions are obtained for the tidal disturbance at 
any point whatever, and at certain particular points, such as points on the 
boundary. The last of the five cases treated is that of a canal lying along a 
parallel of latitude and would appear to furnish a means of checking the same 
case treated by Airy's Canal Theory of Tides. 
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